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Lattice discretisation errors in the Landau gauge condition are examined. An improved gauge fixing algorithm 
in which 0(a 2 ) errors are removed is presented. 0(a 2 ) improvement of the gauge fixing condition displays 
the secondary benefit of reducing the size of higher-order errors. These results emphasise the importance of 
implementing an improved gauge fixing condition. 



1. Introduction 

Gauge fixing in lattice gauge theory simulations 
is crucial for many calculations e.g. the study 
of gauge dependent quantities such as the gluon 
propagator [Q. However, the standard lattice 
Landau gauge condition Q is the same as the con- 
tinuum condition, d^A^ — 0, only to leading 
order in the lattice spacing, a. 

The focus of this talk is to use mean-field- 
improved perturbation theory || to compare dif- 
ferent lattice definitions of the Landau gauge, 
and quantify the sizes of the discretisation errors. 
In particular, we derive a new 0(a 2 ) improved 
Landau-gauge-fixing functional, and a method of 
generalising this to 0{a n ). 

2. Lattice Landau Gauge 

Gauge fixing on the lattice is achieved by max- 
imising a functional whose extremum implies the 
gauge fixing condition. The usual Landau gauge 
fixing functional is [|| 

^[{U}} = £ ±Tr {U°(x) + C/«(.T)t} , (1) 
where U®(x) = G(x)U^(x)G(x + /i) f . By taking 



J(x 

the functional derivative of O), it can be shown 
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that a maximum of that functional implies the 
continuum Landau gauge, with 0(a 2 ) errors. It 
can further be shown that the gauge fixing con- 
dition implies that 



[d^A^x) 



(2) 



where Tii represents 0(a i ) and higher-order 
terms. Naively one might hope that higher-order 
derivatives in the brackets are small, but it will 
be shown that the terms on the R.H.S. of (||) are 
large compared to the numerical accuracy possi- 
ble in gauge fixing algorithms. 

This "one-link" functional can be generalised 
to functionals using "n-link" terms: 



J*[U\ = T i ±Tr{U%(x) + h.c.} 
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where 
Then 
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0(a 2 ) errors can be removed from the gauge 
fixing condition by taking a linear combination of 
the one-link and two-link functionals: 



G 



(5) 



T G - 4 ^G 1 
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where we have included the plaquette-based, 
mean-field (tadpole) improvement parameter, 
§• 

To perform the gauge fixing we adopt a "steep- 
est descents" approach Q]. The gauge transfor- 
mation is G(x) = exp{—ia d^A^ix)}. To 
maximise, for example, J 7 ^ , we use (|4|) to derive 
the gauge transformation 



exp 



Gi(z) 
where 

Ax(^) = [U^i* A) - U°{x) - h.c.} 



(6) 
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Similarly, A2 and Ai mp are obtained from the 
functional derivatives of T% an d Ji mp respec- 
tively. For a given functional, J- G , the gauge 
fixing algorithm proceeds by calculating the rel- 
evant A,, applying the associated gauge trans- 
formation to the gauge field, and iterating until 
the lattice Landau gauge condition is satisfied, to 
within some numerical accuracy. The approach 
to Landau gauge is measured by 



(7) 



A configuration fixed using Ai(x) will satisfy 
(||). It will also satisfy 

,2 



^ 2 E 



A 2 (x) = —2iga 

and similarly, 
Ai mp (x) = -2iga 2 



dMx)-Hi + H 2 



—Hi + Wimp} . (8) 

Since the improved measure has no 0(a 2 ) error of 
its own, (||) provides an estimate of the absolute 
size of these discretisation errors. 



Table 1 

Values of the gauge-fixing measures obtained us- 
ing the improved gluon action on 6 4 lattices at 
three values of the lattice spacing, fixed to Lan- 
dau gauge with the one-link, two-link and im- 
proved functionals respectively. 
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3. Calculations on the Lattice 

We use an 0(a 2 ) tadpole- improved action. For 
the exploration of gauge fixing errors we con- 
sider 6 4 lattices at f3 = 3.92, 4.38, and 5.00, cor- 
responding to lattice spacings of approximately 
0.35, 0.17, and 0.1 fm respectively. 

The configurations are gauge fixed, using Con- 
jugate Gradient Fourier Acceleration Q until 
0\ < 10~ 12 . 0i mp and 02 are then measured, to see 
the size of the residual higher order terms. The 
evolution of the gauge fixing measures is shown 
for one of the lattices in Fig. [l]. This procedure is 
then repeated, fixing with each of the other two 
functionals, and the results are shown in Tabic [lj. 
Results from additional lattices, as well as a more 
detailed discussion, are in ||. 

If we fix a configuration to Landau gauge by 
using the basic, one-link functional, then the im- 
proved measure, #i mp , will consist entirely of the 
discretisation errors. Looking at Table [I], we see 
that at /3 = 4.38, #i mp = 0.059, a substantial de- 
viation from the continuum Landau gauge com- 
pared to the tolerance of the gauge fixing. We 
note that the relationship between the function- 
als in (|^) provides a constraint on the gauge fixing 
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Figure 1. The gauge fixing measures for a 6 lat- 
tice with Wilson action at (3 = 6.0. This lat- 
tice was gauge fixed with Ai, so 9\ drops steadily 
whilst 02 and 9i mp plateau at much higher values. 



measures. For example, when fixing with Ai 



-'Imp 



(- — ) 2 

V 12 J 



12 



1)2 



= - ~ 0.111. 

9 



(9) 



A configuration fixed using Ai mp (a;) will satisfy 
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With the help of (§) we see that 
A 1 (x) = -2iga 2 J2{^Mx) 



Hi — Tii 
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and (H) are identical to within a sign. If the three 
different methods presented all fixed in exactly 
the same way, then the #i mp of a configuration 
fixed with Ai, would be equal to 9\ when the 
configuration is fixed with Ai mp . It is clear from 
the table that they are not, signaling the higher- 
order derivative terms d^A^(x), contained in the 
the TLi, take different values depending on the 
gauge fixing functional used. 

Examining the values in Table [j] reveals that in 
every case 9\ is smaller when we have fixed with 
^imp than flimp under the T\ . This suggests that 
the additional long range information used by the 
improved functional is producing a gauge fixed 
configuration with smaller, higher-order deriva- 
tives; a secondary effect of improvement. 



Equally, one can compare the value of 02 when 
fixed using T\, and 9\ when fixed using the JT 2 . 
In this case, their differences are rather large and 
are once again attributed to differences in the size 
of higher-order derivatives of the gauge field. Ti 
is coarser, knows little about short range fluctu- 
ations, and fails to constrain higher-order deriva- 
tives. Similar conclusions are drawn from a com- 
parison of 92 fixed with the Ji mp and #i mp fixed 
with J~2- 

We also find that in terms of the absolute er- 
rors, the Wilson action at (3 — 6.0 is comparable 
to the improved lattice at /? = 4.38, where the 
lattice spacing is three times larger. 

4. Conclusions 

We have fixed gluon field configurations to Lan- 
dau gauge by three different functionals: one-link 
and two-link functionals, both with 0(a 2 ) errors, 
and an improved functional, with C(a 4 ) errors. 
Using these functionals we have devised a method 
for estimating the discretisation errors involved. 
Lattice Landau gauge, in its standard implemen- 
tation, deviates from its continuum counterpart 
by one part in 20, despite fixing the Lattice gauge 
condition to one part in 10 12 . Our results indicate 
that order 0(a 2 ) improvement of the gauge fix- 
ing condition improves comparison with the con- 
tinuum Landau gauge through: 1) the elimina- 
tion of 0(a 2 ) errors and 2) reducing the size of 
higher-order errors. 
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